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ABSTRACT 

The Auger observatory has observed a possible correlation between Ultrahigh Energy Cosmic Rays 



Qh (UHECRs) above 57 EeV and nearby candidate Active Galactic Nuclei (AGN) from the Veron-Cetty 

Veron catalog (VCV). In this paper we report on Chandra X-ray observations of 10 unconfirmed VCV 
"^2 ' AGN-candidatcs and luminous IR galaxies correlating with the first set of Auger UHECRs, to deter- 

mine whether or not they have active nuclei. The X-ray data, when combined with optical luminosities, 
show that in fact none of the 10 galaxies have a significant AGN component; if there is any nuclear 
| activity at all, it is weak rather than obscured. This reduces the number of UHECRs in the original 

, Auger dataset possibly correlating with AGNs from 20 of 27 down to 14 of 27. We also used Chandra 

(N ■ 

to measure the X-ray luminosity of ESO 139-G12, an AGN which correlates with 2 of the Auger 
Q\ • UHECRs, to obtain the first estimate of its bolometric luminosity; this completes the determination 

of the bolometric luminosities of all correlating AGNs. Taking our results into account, only one of 
the 27 UHECRs in the original Auger data-release is correlated on a few-degree angular scale with an 
identified AGN that is powerful enough in its steady-state to accelerate protons to the observed en- 
ergies, according to conventional acceleration mechanisms. Intriguingly, w 60% of the UHECRs with 
\b\ > 10° do correlate with genuine AGNs, but these are too weak to meet the acceleration criterion 
for protons; this may be an indication that AGNs experience transient high-luminosity states which 
can accelerate UHECRs. To determine the source(s) and composition of UHECRs through statistical 
correlation studies requires reliable, complete and uniform catalogs of identified AGNs; our a poste- 
riori inspection of ambiguous source candidates underscores the inadequacies of the VCV catalog in 
this respect. 

1. INTRODUCTION of high energy events and because the deflections of the 

Identifying the sources of ultra-high energy cosmic rays cosmic v ^ s as the y travel throu S h intergalactic magnetic 
(UHECRs) is one of the major outstanding goals in as- fields mean that thc cosmic Y ^ arrival directions do not 
trophysics. Progress has proven difficult due to the rarity P oint back to their ori S ins - The ver y most energetic cos- 

mic rays (E > 6 x 10 19 eV) have attracted attention as a 



promising path forward. The energy loss due to the GZK 
mechanism means that CRs with such energies typically 
have traveled about 100 Mpc, significantly limiting the 
possible sources for such high energy particles. Further- 
more, at these energies magnetic deflections of protons 
may be small enough to allow the identification of the 
progenitor type based on statistical associations. 

By 2007, the Pierre Auger Collaboration (2007, 
2008) had compiled enough events to begin draw- 
ing statistical correlations between UHECR events 
and possible sources, and reported a strong cor- 
relation betw een UHECRs and the nearby galaxies 



expla ined by nearby galaxy clustering alone (Za w et al 



listed in the IVeron-Cettv and Veronl ([20061 ) Catalog 
of Quasars and Active Galactic N uclei (12 th Ed.) 
( The Pierre Auger Collaborationll2007l ). The scan proce- 



dure that was used steps through UHECR energy thresh- 
old, maximum angular separation, and maximum VCV 
galaxy redshift to find the values giving the lowest chance 
probability compared to an isotropic distribution, then 
evaluates the likelihood of such a correlation occuring 
by chance by performing the same analysis on many 
isotropic datasets. In the following, the term "corre- 
lated" referring to a galaxy with respect to a UHECR 
simply means that the given galaxy falls within the an- 
gular and redshift limits returned by applying the scan 
method with that galaxy catalog. 

Of the 27 cosmic rays with energies above 57 xlO 18 eV 
that Auger detected before Aug. 31, 2007, twenty corre- 
late within 3.2° with VCV galaxies having z < 0.018 , 



about 75 Mpc ([The Pierre Auger Collaboration 



2008) 



Since more than one galaxy can be correlated with a 
UHECR and vice versa, there are 21 VCV galaxies cor- 
related with these 20 UHECRs. Galaxy catalogs such as 
VCV are incomplete in the Galactic Plane, so a better 
comparison is obtained by restricting to \b\ > 10° where 
the VCV catalog is more complete. With this restriction 
there are 22 UHECRs of which 19 UHECRs correlate 



with 20 galaxies (|Zaw et al.ll2009l) . This is a much higher 
correlation than would be expected by chance from an 
isotropic source distribution, and higher than can be 



20111 ). More recent Auger data continue to show a signif- 



icant, albeit less stron g, correlation with VCV galaxies 
(|Abraham et al.ll2009l ). 

An important question is whether the observed cor- 
relation with VCV galaxies implies that AGNs are the 
sources of some or all UHECRs. The VCV catalogue 
is a list of AGN candidates, so first it must be estab- 
lished whether the correlating galaxies (i.e., the VCV 
galaxies within 3.2° of a UHECR) are actually AGNs 



or not. 



Zaw et al 



(2009) looked at existing observa- 
tions of the galaxies in VCV that were correlated with 
the first set of UHECRs and found that only 14 of the 
21 galaxies show unambiguous evidence of AGN activ- 
ity. Three show no signs of AGN activity, while the 
other four have ambiguous optical spectra. A widely 
used technique for optical identification of AGNs, so- 
called BPT line ratios, com pares the relative str engths 



of diagnostic spectral lines ((Baldwin et al 



1981 ). The 



BPT line ratios of the 4 ambiguous galaxies f all within 
the ranges classified bv lKauffmann et al] ([20031 ) as AGNs 



bu t they do not f a ll wit hin the line-ratio ranges adopted 
by iKewlev et al] (|200l[ ) as indicative of an active nu- 
cleus. AGN activity can often be obscured in the optical 
bands by dust obscuration, and the Kewley test excludes 
some known AGNs. Additionally, as many as half of 
AGNs that are selected based on radio or X-ray proper- 
ties would not be identified by looking only at their BPT 



line ratios ( Reviglio and Helfand 



2006). Thus, determin- 



ing whether these 4 ambiguous VCV galaxies have ac- 
tive nuclei requires observations outside of optical wave- 
lengths. 

The VCV galaxies are not the only population of 
nearby galaxies foun d to co rrelate with the Auger UHE- 
CRs. 



Berlind et al 



(|2010f ) find an anomalously large 
correlation between the Auger galaxies and nearby Lu- 
minous IR galaxies in the PSCz catalog. Restricting 
to |6| > 10° where PSCz is complete, 13 galaxies of 
Lib. > 10 10 ' 5 £ S im correlate within 2.1° with one or more 
of the 22 Auger UHECR events. Some LIRGs contain 



an active nucleus with dust absorbing the AGN radiation 
and re-emitting it thermally, giving rise to large IR lumi- 
nosities. This raises the question of whet her the LIRGs 
found to correlate may actually be AGNs 



Berlind et al 



(|2010f) showed that 6 of the 13 correlating IR galaxies 
are also in VCV and 5 of these do in fact host AGNtQ- 
The other 7 correlating IR galaxies lacked the observa- 
tions needed to differentiate between star-formation and 
obscured nuclear activity. 

A key distinguishing feature of active galaxies is that 
accretion-driven radiation produces a nearly flat spec- 
trum from the IR to the X-ray, known as the broad- 
band continuum. Normal and starburst galaxies, on the 
other hand, have broadened blackbody spectra peaked 
sharply in the UV/optical. This means the X-ray to op- 
tical flux ratios are considerably harder for AGNs than 
for normal/starburst galaxies. If the active nucleus is ob- 
scured, the X-ray to optical ratio becomes even harder 
since dust ab sorbs UV and optical photons more readily 



than X-rays (|Comastri et al 



20031) . This makes X-ray 



observations, particularly when combined with observa- 
tions in the near-IR, optical or UV, a powerful tool for 
identifying obscured AGNs. 

In this paper we use X-ray observations to determine 
whether the 4 ambiguous VCV galaxies, and the 7 in- 
determinate PSCz galaxies, have active nuclei. We ob- 
served ten of these possible UHECR source galaxies using 
the Chandra X-ray satellite, while for IC 5169 we used 
data from a recent XMM-Newton observation. 

Another interesting question is whether AGNs that 
correlate with UHECRs have any characteristic features 
which distinguish them from the AGNs that do not seem 
to correlate with UHECRs. A particularly relevant prop- 
erty is the luminosity of the AGN, as this helps to con- 
strain the possible cosmic r ay acceleration mechanisms 



(Farrar and Gruzinov 



2008). Previously, reliable lumi- 



nosities have been established for all but one of the cor- 



relating AGNs (|Zaw et all 120091) . (jBerlind et al 



2010). 



1 The one which is not a confirmed AGN falls within the VCV 
sample to be tested. 



We observed the remaining AGN, ESO 139-G12, in or- 
der to obtain the first robust estimate of its bolometric 
luminosity. 

2. OBSERVATIONS AND DATA REDUCTION 

Each target was observed on the Chandra ACIS-I de- 
tector for 5 kiloseconds between January and August 
2009. Data analysis was performed using the ciao data 
analysis software. All of the targets corresponded to an 
X-ray source. (Table [TJ. 

For the Chandra observations, we used two different 
methods to estimate the X-ray energy flux. The direct 
method sums the energy received per unit time in the 
band of interest (2-10 keV) in a 2 arcsec window around 
the source using the eff2evt tool which takes into ac- 
count both the quantum efficiency and the effective area 
corrections of the satellite. The flux is typically domi- 
nated by higher energy and less frequent events, making 
this technique vulnerable to Poisson fluctuations when 
the count rates are small. The second technique con- 
siders only the total counts received from the compact 
nuclear source, ignoring the energies of the photons, and 
assumes that the spectrum follows a power law with a 
spectral index of 2. Then webPIMMS can estimate the 
corresponding flux, including Galactic extinction from 
webPIMMS and assuming no intrinsic extinction. For 
the galaxies which are diffuse sources, we report only the 
direct measurement of the flux within 2 arcsec of the 
galactic center, as an upper limit on the possible nuclear 
emission. 

We also analyzed the XMM-Newton observation of IC 
5179. Using the source imaging PPS processing provided 
by XMM, we took the count rate on the EPIC pn cam- 
era (0.5-12 keV) and used webPIMMS to estimate the 
corresponding 2-10 keV flux of the source. 

3. ANALYSIS AND RESULTS 

X-ray to optical flux ratios provide a useful way 
to classify galaxies according to their nuclear activity. 



Barger et al 



(2002) show that for AGNs, the ratio of 



the 2-10 keV flux to the R-band magnitude is typi- 
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ESO 565-G006 
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5.5 E -14 


2.4 E -14 


12.7 1 


-2.7 


3.2 E 40 


NGC 7648 
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3.0 E -15 
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-4.4 
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2MASX Jl 754-60 


10262 
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2.7 E -15 








1.6 E 39 


IC 5179 


XMM-Newton 


2685 






9.1 E -14 


11.4 


-2.98 


2.5 E 40 


ESO 139-G12 


10263 


1070 


666 


4.7 E -12 c 




12.9 1 ' 2 


-0.2 


2.9 E 43 



TABLE 1 
Table listing source properties 



1 20051 . (3l lde Vaucouleurs and Loneol 11198811 



References. — (11 lLauberts and Valentiinl lll989ll . f2l iDovle et al. 
Note. — Fluxes are for 2-10 keV and are in cgs units (erg cm" 2 s" 1 ); only the direct flux determination method is used for diffuse 
sources. The 2-10 keV luminosity, L2— i0i is in units of erg s — 1 and is derived from the direct flux measurement; it is the total X-ray 
luminosity in the central region and thus an upper limit on X-ray luminosity of a possible AGN. The first four rows are the correlating 
galaxies found in VCV. The following six are those from PSCz. The next row was observed by XMM-Newton: the counts are from 0.5 - 
12 keV and the observation lasted 25000 sees. The last row is a known AGN whose X-ray luminosity had not previously been measured, 
needed to assess its ability to accelerate UHECRs. 

(a) No counts above 2 keV. (ft) Johnson magnitude rather than Cousins. This does not affect the result, (c) Derived from fit, see Figure 1. 



cally 1 > ^o(7(/ I //fl]t]> —1, while obscured AGNs have 
l°9(fx/ fa) > 1- Normal and starburst galaxies have 
l°9(fx/ Ib) < ~ 2, making this a robust way of discrimi- 
nating AGNs, even obscured ones, from normal and star- 



burst galaxies (|Donlev et al 



2005) 



The results are summarized in Table [TJ The X-ray to 
optical ratios are far below what is expected from AGNs 
(log(f x / fn) > —1) for ten of the ambiguous or indeter- 
minate galaxies: nine we observed with Chandra, and IC 
5179 for which we used an XMM-Newton observation. 
This holds even when the larger of the two X-ray flux 
determinations is used to calculate the X-ray to optical 
flux ratios. Since our measure of nuclear activity depends 
on logio(f x ), h is relatively insensitive to the uncertain- 
ties inherent in the determination of the flux at these 
low count rates. We conclude that none of these galaxies 



f l°g(fx/.fn) = log\o{fx) + 5.5 - 
keV flux 



.4 * Rmag where f x is the 2-10 



have active nuclei. 

The last of the indeterminate galaxies (2MASX Jl 
754- 60) lacks the optical observations needed for the 
logifx/ fa) test, but its X-ray emission is diffuse and 
weak, and its near-IR properties are similar to the other 
galaxies. Therefore we conclude that it is unlikely to 
host an AGN, certainly not one expected to be capable 
of UHECR acceleration. 

Obscured AGN in X-ray surveys are often identified by 
the hardness ratio of the source. The hardness ratio is 
defined as H-S/H+S where H is counts from 2-8 keV and 
S is counts from 0.5-2 keV. Typically, nearby AGN with 
column densities greater than Nh ~ 10 22 have positive 
hardness ratios, while unobscured sources have negative 



hardness ratios ( Fiore et al 



20001 ) . Table Q] gives counts 



from 0.5-10 keV and 2-10 keV. These can be used to es- 
timate the hardness ratio since the small effective area 
of Chandra at high energies meant that no photons with 
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Fig. 1.— Data and fit to a power law (r = 0.72+°-[5g, 90% 
confidence) with Nh < 8.2 X 10 21 cm -2 at 90% confidence level. 
This is our best fit model of the data, with a reduced \ 2 of 1.7. 
This fit corresponds to a 2-10 keV flux of 4.7±°-f x !0~ 12 erg cm -2 
s — 1 (error denotes 90% confidence interval). 

energies above 8 keV were detected for any of the sources. 
Although the small number of counts means large uncer- 
tainties in the hardness ratios, they show no evidence of 
obscuration. 

4. LUMINOSITY OF ESO 139-G12 

The AGN ESO 139-G12 is correlated with two UHE- 
CRs, making it a particularly interesting source candi- 
date. This source is bright enough to determine its flux 
by fitting the spectrum using xspec. The fit, assuming 
an absorbed power law, returned a hard photon index 
of 0.72^0 09 (error denotes 90% confidence interval) and 
very low absorption, Nh < 8.2 x 10 21 cm~ 2 at 90% con- 
fidence level (Figure 1). From this fit we extracted a 
2-10 keV flux of 4.7±?;? x 10~ 12 erg cm" 2 s" 1 (error de- 
notes 90% confidence interval). When combined with its 
redshift, z = 0.017 (~ 74 Mpc assuming a standard cos- 
mology with H Q — 70), we find its 2-10 keV luminosity: 
£2-10 ~ 3.liQy x 10 42 ergs -2 . The conversion factor 
relating £2-10 to the bolometric luminosity depends on 
the r elative activity of the AGN (jVasudevan and Fabian 
20071) . In the absence of a black hole mass measure- 
ment for this galaxy, we follow 
Lboi = 20 L2 -10 as is typical for an AGN with Eddington 
ratio less than 0.1. This gives Lboi 



Zawetal 



(120091 ) and take 



of ESO 139- G '12. 

Although ESO 139- G 12 is the only galaxy we observed 
which appears to have nuclear activity, it is always possi- 
ble that the host galaxy simply dominates the AGN ener- 
getically, making the AGN impossible to detect. There- 
fore the 2-10 keV luminosity is given for the other galax- 
ies as well. With an appropriate conversion factor, this 
can be used to get an upper limit on possible weak nu- 
clear activity in the galaxy which may be of interest in 
other studies of UHECR sources. 

5. CONCLUSIONS 

These Chandra observations complete the task of 
determining which galaxies found to correlate with 
UHECRs bybhe Pierre Auger Collaboration! (|2007l ) and 



Berlind et al 



6.2+? 8 x 10 4 



erg 



the first measurement of the bolometric luminosity 



(2010) have active nuclei, and of determin- 
ing the bolometric luminosities of the AGNs correlating 
with UHECRs. The X-ray fluxes of the ten unclassified 
correlating galaxies we studied all fall within the range 
typical of normal and starburst galaxies of the same opti- 
cal magnitude, hence we find no evidence of AGN activity 
in any of themj 2 ] 

Thus only 14 out of 27 UHECRs (13 out of 22 UHECRs 
with |6| > 10°) correlate within 3.2° with an actual AGN, 
using the Auger scan method to correlate UHECRs with 
VCV galaxies but discarding candidate sources which are 
not in fact AGNs. Of the 13 highly luminous IR galaxies 
{Lm, > 10 115 i©) found to correlate within 2.1° with 
one (or more) of the 22 UHECRs with \b\ > 10°, we find 
that only five are also AGNsi Five of the 27 UHECRs, 
and one of the 22 UHECRs with |6| > 10°, have neither 
an AGN nor a LIRG within 3.2° or 2.1° respectively. 
These uncorrelated UHECRs do not preferentially have 
low energies, for which sources are typically more distant, 
so the lack of correlation cannot be attributed entirely 
to incompleteness of the VCV or PSCz catalogs at larger 

2 Of course, it is impossible to rule out the possibility of very low 
luminosity nuclear activity which is energetically dominated by the 
host galaxy, but such weakly active AGN are not well-motivated 

candidates for UHECR acceleration. 

3 We could not run this correlation for the full 27 UHECRs 

because IRAS does not observe within the Galactic plane. 
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redshifts. 

The degree of correlation with VCV galaxies was 
found to be lower in the second Auger data release 



( Abraham et al 



20091) . It has also been shown recently 



that the scan method itself is subject to large sample-to- 
sample variance and that it systematically tends to over- 
estimate the degree of correlation (Farrar et al., 2011, in 
preparation). These facts may suggest that other can- 
didate sources not associated preferentially with AGNs 
may be responsible for some, most, or all UHECRs. But 
other factors can contribute to a lack of correlation be- 
tween UHECRs and AGNs as well, even if active galaxies 
are an important source of UHECRs: incompleteness of 
source catalogs (especially within the Galactic plane), 
magnetic deflection that obliterates angular correlation 
between the CRs and their true source (particularly for 
UHECRs with charge Z > 1), or transient production 
of UHECRs via exceptional, po werful flares flares in 
very weakly or non-active galaxies (jFarrar and Gruzinov 
200J). 

Another possibility raised by our results is that the 
VCV catalog may be so badly contaminated with non- 
AGNs, that the correlation anal y ses ap plied to the data 
up to now, c.f.. 



Abraham et al 



([20091 ) . may be thrown 



off. The "spot check" of VCV provided by this work 



combined with (Zawetal 



20091 ) is the most thorough to 



date, by providing definitive identifications of AGN can- 
didates; if it gives a representative picture of the full VCV 
catalog, a third of the galaxies in the VCV catalog are 
not true AGNs. Earlier spot-checks using (incomplete) 
data available in the literature, but considering a larger 
angular region around UHECR s and/or later UHECR 

|201lL find that 1/4 - 1/3 of 



Zaw et al 



catalogs, e.g. 

the AGN candidates in VCV are not AGNs. The impact 
of a contamination of this magnitude on statistical cor- 
relation analyses such as Auger's scan method needs to 
be examined with simulations - it may bias the results 
in unpredictable ways. 
The one established AGN which we observed in order 



to determine its bolometric luminosity, ESO 139-G12, 
proves to have Lboi comparable to m ost of the other 



20091) . Knowledge of 



13 correlated AGNs ([Zaw et al. 
Lboi and the energies of the correlated UHECRS (89 and 
59 EeV) allows us to evaluate Ab oi = 10~ 4 5 £boi ffi opi 



Zaw et al 



(2009J) to 



the figure-of-merit introduced by 
quantify the ability of an AGN to accelerate a proton 
to the energy of the correlated UHECR. A value of 
Aboi 1 satisfies the acce l eratio n criterion for protons 
(c.f., iFarrar and Gruzinovl (|2009l )). With A bo i = 0.1, 
ESO 139-G12 is weak for standard UHECR acceleration 
mechanisms. 

6. SUMMARY 

Our observations and analysis using the Chandra X-ray 
satellite and other data establish that one-third of the 21 
galaxies in the Veron-Cetty Veron catalog of AGN can - 
didates found bv lThe Pierre Auger Collaboration! (|2007l ) 
to correlate with UHECR arrival directions, do not in 
fact have active nuclei. Combining this with our mea- 
surement of the X-ray luminosity of ESO 139-G12, an 
AGN correlating with 2 UHECRs, implies that only one 
of the 27 UHECRs in that first Auger data-release corre- 
lates with an AGN that is powerful enough in its steady- 
state to accelerate protons to the observed energies, ac- 
cording to conventional acceleration mechanisms. How- 
ever 60% of the UHECRs with \b\ > 10° do correlate with 
genuine but weak AGNs, consistent with the possibility 
that AGNs may have a flaring state th at can acceler- 
ate UHECRs (^rrarandGruzin3[2009). These results, 



along with the second Auger data release showing a lowe r 



correlation even with weak AGNs (jAbraham et al.ll2009l) , 
underscore the fact that the source(s) of the highest en- 
ergy cosmic rays, as well as their composition, remains 
an open question. 
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and NSF-PHY-0900631. GRF and IZ acknowledge their 
membership in the Pierre Auger Collaboration. 
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